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ABSTRACT 


The  clinical  comparison  of  two  treatments  poses  an  ethical 
problem  for  the  experimenter,  who  wishes  to  minimize  the  number  of 
assignments  of  the  poorer  treatment  while  maintaining  statistical 
significance.  Sequential  methods  using  data-dependent  treatment 
allocation  offer  one  possible  route  out  of  this  dilemma.  We  give 
a  general  formulation  of  two-population  sequential  testing,  and 
allocation  rules  for  the  normal  and  exponential  distributions. 

Monte  Carlo  trials  indicate  that  a  new  rule  for  the  exponential  case 
is  superior  to  previously  known  rules. 
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Chapter  I  Sequential  Testing  and  the  Clinical  Trial 


A.  Introduction 

One  of  the  most  fertile  sources  of  statistical  problems  has 
long  been  the  controlled  clinical  trial.  We  do  not  propose  to  give  a 
detailed  history  of  this  field  here:  this  would  be  a  major  project 
in  its  own  right.  Rather,  we  will  concentrate  our  attention  upon  the 
use  of  sequential  methods  in  clinical  settings,  and  some  of  the  problems, 
both  practical  and  statistical,  arising  therefrom.  In  particular,  we 
shall  consider  some  methods  of  sample-size  reduction  in  two-population 
hypothesis  testing.  Sequential  estimation  is  not  discussed. 

Sequential  analysis  can  be  defined  as  that  area  of  statistics 
which  assumes  the  sample  is  drawn  unit  by  unit  in  time.  After  the 
selection  of  each  unit  a  decision  is  made  whether  or  not  to  continue 
sampling,  and  if  not,  to  make  one  of  a  number  of  decisions.  In  its 
early  years,  sequential  analysis  was  largely  the  work  of  Abraham  Wald, 
whose  book  (19),  remains  the  classic  in  the  field.  The  use  of 
sequential  methods  is  suggested  by  many  experimental  situations,  par¬ 
ticularly  those  in  which  complete  control  of  the  sample  is  not  possible. 
This  is  manifestly  the  case  in  many  medical  situations.  For  example, 
one  cannot  usually  select,  say,  fifty  candidates  for  a  new  cancer 
treatment:  the  patients  must  come  to  the  doctor  first,  providing  a 

natural  basis  for  sequential  analysis. 

Medical  experimentation  has  a  dimension  which  it  shares  with 
few  other  fields  -  an  ethical  one.  Weinstein  (20)  discusses  this 
problem,  and  gives  a  short  bibliography.  Statistical  requirements  of 
sample  size  often  conflict  with  medical  goals.  A  sequential  design 
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for  an  experiment  is  immediately  attractive  to  the  medical  researcher 
since  it  offers  the  possibility  of  curtailing  administration  of  an 

l 

inferior  treatment  much  sooner  than  would  be  possible  with  a  fixed- 
size  design.  For  this  reason,  medical  researchers  were  early  attracted 
by  Wald's  methods.  The  early  developments  along  these  lines  are  well 
documented  in  Armitage  (3),  (4)s  and  Anscombe  (2). 

The  sequential  probability  ratio  test  (SPRT)  is  the  centre¬ 
piece  of  Wald's  work.  However,  most  clinical  researchers  have  rejected 
it  as  unsuitable  for  medical  purposes.  To  better  appreciate  this,  we 
now  give  a  brief  discussion  of  this  test. 

B .  The  Sequential  Probability  Ratio  Test 

Let  (ft,  B,  P)  be  a  probability  space  on  which  is  defined  a 

OO 

sequence  (X  )  of  independent,  identically  distributed  random  vari- 
n  n=l 

ables  which  have  density  p(x;0)  with  respect  to  P  for  0  £  H  .  To 
sequentially  test  the  two  hypotheses 

Ho  :  0  =  0O 

H1  =  0=  ®1  V- 

two  constants  a  and  b  are  chosen  such  that  0  <  a  <  1  <  b  <  °°.  is 

ttl 

observed,  followed  by  X^,  X^,  and  so  on.  After  the  n  observation 
(n  =  1,  2,  ...)  the  probability  ratio  L  is  calculated. 

jr  P  (Xi;  0^ 

L  =  - 

n 

jr  P  (Xi;  0Q) 


i  =  1 
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The  following  decision  rules  are  used: 


1. 

If 

L  '  b , 
n 

terminate 

sampling 

and 

accept  H^. 

2. 

If 

L  <  a, 
n 

terminate 

sampling 

and 

accept  H  . 

0 

3. 

If 

a  —  Ln  — 

b ,  observe  X  , _  . 

n+1 

The  constants  a  and  b  are  chosen  such  that 
P  (Accepting  |  H  )  £_  a  ,  and 
P  (Accepting  Hq  |  H^)  3. 

If  is  designated  as  the  null  hypothesis,  then  a  is  the  probability 
of  Type  I  error,  3  the  probability  of  Type  II  error.  As  in  the  fixed- 
sample  tests,  a  and  3  are  pre-assigned. 

1 .  Properties  of  the  SPRT 

We  now  state,  without  proof,  the  four  most  important  properties 
of  the  SPRT.  For  proofs,  we  refer  the  reader  to  Wald. 

1.  The  test  terminates  with  probability  one.  That  is,  a  decision  is 
eventually  made. 

2.  b  <  1  -3  and  a  >  3 

a  1  ~  a 

3.  Let  N  =  min  (n  :  £  (a,b));  that  is,  N  is  the  size  of  the  sample 

at  termination.  Then  for  all  m  <  °°,  P  (N  >_  m)  >0. 

4.  Among  all  tests  of  H  against  with  the  same  error  probabilities, 

the  SPRT  minimizes 

E  (N  |  0  )  and  E  (N  I  0J  . 
o  1 

The  fourth  property  is  the  test’s  most  medically  attractive 
point.  It  is  the  third  which  renders  the  SPRT  unsuitable  for  clinical 
use.  No  doctor  is  likely  to  countenance  the  possibility  of  systematically 
giving  a  poor  treatment  to  an  indefinitely  large  number  of  patients. 
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For  this  reason  much  attention  has  been  given  to  methods  of  truncating 

or  restricting  sequential  tests  in  order  to  reduce  sample  sizes.  Notable 

/ 

advances  in  this  domain  were  made  by  Armitage  (3)  and  Anderson  (1). 

C .  Sample-Size  Problem  in  the  Two-Sample  Case 

In  the  foregoing  discussion  we  have  assumed  that  there  is 
only  one  treatment  or  population  in  question.  When  two  or  more  are  to 
be  compared,  the  desiderata  for  an  experimental  protocol  are  somewhat 
different.  In  the  one-population  case,  one  typically  attempts  to 
determine  if  a  treatment  is  acceptable  according  to  some  established 
criterion.  With  several  competing  treatments,  the  experimenter  more 
often  wishes  to  select  the  "best",  or  to  rank  them.  In  this  setting, 
the  overall  sample  size  becomes  less  crucial.  Instead,  it  is  desired 
to  minimize  the  number  of  times  the  poorer  treatments  are  administered. 

In  this  thesis,  we  are  concerned  with  the  two-sample  case.  For 
a  comprehensive  treatment  of  k-sample  identification  and  ranking  methods, 
see  Bechhofer,  Kiefer  and  Sobel  (.6)  ,  who  consider  only  vector-at-a-time 
and  cyclic  sampling,  ignoring  the  sample  size  problem  we  are  considering. 
The  test  we  will  use  is  a  generalization  of  the  SPRT:  its  discussion 
is  deferred  to  the  next  chapter. 

D.  The  Bechhofer-Kiefer-Sobel  Procedure  P  * 

15 

The  basic  ranking  procedure  developed  by  Bechhofer,  Kiefer  and 
Sobel  for  k  populations  has  as  its  primary  goal  the  selection  of  the 
set  of  the  t  "best",  where  t  could  be  anything  from  1  to  k-1.  We  shall 
assume  t  =  1.  Before  any  observations  are  taken  it  is  assumed  that 
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each  treatment  is  equally  likely  to  be  chosen.  After  each  vector  of 
observations  has  been  made,  a  posterior  probabilities  of  correct 

/ 

selection  are  calculated.  The  first  time  one  of  these  probabilities 
exceeds  the  desired  probability  of  correct  selection,  the  procedure 
terminates  and  the  treatment  corresponding  to  this  maximum  probability 
is  selected. 

E.  Types  of  Response  and  Related  Probability  Models 

Medical  testing  situations  fall  roughly  into  three  broad 
categories,  each  characterized  by  what  we  shall  refer  to  as  a  "response 
type"  and  an  associated  archetypical  distribution.  Each  type  has 
attracted  considerable  attention  as  a  setting  for  sequential  trials. 

1.  General  Continuous  Response 

In  many  cases,  patient  response  to  a  treatment  is  such  that  it 
allows  the  assumption  of  continuity  of  the  response  variable.  Two 
examples  are: 

a)  Blood  pressure  change  following  administration  of  a 
hypotensive  agent,  and 

b)  Weight  change  following  a  fixed  period  of  supervised  diet. 
The  Normal  distribution  enters  naturally  at  this  point,  if  for 

no  other  reason  than  the  Central  Limit  Theorem.  Many  authors  have 
considered  sequential  tests  for  normal  parameters,  and  discussion  of 
this  area  is  beyond  the  scope  of  this  thesis.  We  will,  however, 
examine  some  recent  results  pertaining  to  the  sample-size  problem. 

2 .  Discrete  Response 

In  the  case  of  responses  falling  naturally  into  a  number  of 
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categories,  a  discrete  probability  model  is  appropriate.  The  most 
important  case  by  far  is  that  of  dichotomous  response,  with  the 

I 

binomial  distribution  the  related  model.  Examples  of  dichotomous 
response  occur  wherever  treatment  results  in  either  success  or  failure, 
clearly  defined. 

Sequential  procedures  for  selecting  the  best  of  a  set  of  bi¬ 
nomial  populations  have  had  much  attention  in  the  last  decade.  One 
sampling  rule  that  effectively  reduces  inferior  treatment  sample  sizes 
is  the  so-called  MPlay-the-Winner "  rule  (PW) .  Briefly  stated,  PW 
sampling  in  the  two-sample  case  follows  each  success  by  a  repeat  of 
the  treatment,  and  each  failure  by  a  change  of  treatments.  This  rule 
was  pioneered  by  Zelen  (21),  and  Sobel  and  Weiss  (17).  Hoel  and  Sobel 
(13)  compared  PW  sampling  with  several  other  rules,  with  especially 
favorable  results  when  both  success  probabilities  are  high. 

3 .  Survival  Response 

In  terminal  illnesses  such  as  cancer  or  heart  disease,  the 
desired  response  to  treatment  is  prolongation  of  life,  usually  within 
an  observable  time-span,  such  as  three  to  five  years.  Of  the  many  life 
distributions,  the  simplest  and  most  important  is  the  exponential.  There 
is  some  empirical  justification  for  the  use  of  this  distribution  as  a 
model  for  survival  in  terminal  cancer  (  Zelen  (22) ),  although  the 
presence  of  non-terminal  cases  in  the  sample  may  invalidate  the  model 
to  some  degree  (Boag,  discussion  of  Armitage  (3)). 

The  feature  that  distinguishes  survival  response  from  the 
general  continuous  case  is  the  natural  censoring  of  the  data.  Full 
information  on  any  patient  consists  of  the  time-span  from  treatment 
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until  death.  Throughout  the  course  of  an  experiment,  one  must  "make 
do"  with  partial  information,  in  the  form  of  survival  times.  In  this 

I 

thesis,  we  are  primarily  concerned  with  the  sample-size  problem  in  a 
survival  response  situation,  with  the  assumption  of  exponentially 
distributed  survival  times. 

F .  Summary  of  the  Thesis 

In  the  next  chapter,  we  discuss  a  general  formulation  for 
sequential  experiments  in  the  two-sample  case.  Included  therein  is  a 
discussion  of  the  sequential  likelihood  ratio  test,  and  a  description 
of  treatment-allocation  rules. 

Chapter  Three  develops  a  two-sample  test  for  comparing  two 
normal  means,  with  attention  given  to  results  by  Robbins  and  Siegmund 
(16)  and  Flehinger  and  Louis  (9) .  Chapter  Four  discusses  the  cor¬ 
responding  test  for  the  scale  parameters  of  two  exponential  distributions. 
Several  allocation  rules  are  presented,  among  them  one  which  appears 
to  be  superior  to  the  Ry  rules  of  Flehinger  and  Louis  (10)  .  The  com¬ 
parisons  were  made  by  means  of  Monte  Carlo  trials,  the  results  of  which 
are  contained  in  the  tables  ending  chapter  Four. 

The  appendix  contains  the  computer  programs  used  in  the  Monte 
Carlo  trials,  with  full  documentation. 
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Chapter  II  Two-Population  Sequential  Testing 


In  this  chapter,  we  present  a  general  formulation  for 
sequential  testing  when  there  are  two  populations,  and  when  data- 
dependent  treatment  allocation  is  used.  In  addition,  we  will  define 
criteria  according  to  which  allocation  rules  can  be  composed. 

In  a  sequential  clinical  trial,  the  protocol  consists  of  three 
rules  as  follows: 

(i)  Admission  rule:  This  is  the  identification  of  patients  eligible 
for  admission  to  the  experiment,  which  is  essentially  a  medical 
problem.  In  the  sequel,  we  shall  not  deal  with  the  rule, 
except  to  assume  that  admission  is  not  dependent  upon  the  course 
of  the  experiment . 

(ii)  Allocation  rule:  Once  a  patient  has  been  admitted  to  the 

experiment,  he  must  be  assigned  one  or  the  other  treatment. 

Rules  by  which  this  decision  are  made  are  the  primary  topic  of 
this  thesis. 

(iii)  Termination  rule:  At  some  point  the  admissions  must  be  ceased 
and  a  decision  made  on  the  relative  merits  of  the  treatments. 
Almost  all  termination  rules  are  related  to  the  SPRT,  discussed 
in  Chapter  I. 

A.  The  Hypothesis  to  be  Tested 

Suppose  that  we  wish  to  compare  treatments  "l"  and  "2",  and 
that  patient  responses  have  densities  f  and  f^*  These  densities  are 
assumed  to  belong  to  a  parametric  family 

F  =  {  f  (x  I  0) ,  0  e  T  ,  x  c-  I  c  R  } 
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We  assume  T  c  .  Let 

D={A:  0  £  T.  ,0+Aet,  0  -  A  £  T  }  . 

i 

Let  f  £  F,  0  e  T  ,  and  A  /  2  £  D.  We  define 

(x)  =  f  (x  |  0  +  A  /  2) ,  and  (2.1) 

f2  (x)  =  f  (x  |  0  -  A  /  2).  (2.2) 

In  a  three-decision  problem,  we  wish  to  choose  one  of  the  following 
hypotheses : 

Hq  :  A  =0;  i.e.,  the  treatments  are  equal, 

:  A  =  A*;  i.e.,  treatment  1  is  better, 

H2  :  A  =  -  A  *;  i.e.,  treatment  2  is  better.  (2.3) 

Here  A  *  >  0  is  a  constant  representing  a  medically  significant 
difference. 

Sometimes  the  hypothesis  H  may  not  be  of  interest,  and  we  wish 

o 

to  choose  one  of  or  only.  We  shall  refer  to  this  latter  formula¬ 
tion  as  the  "two-decision"  problem. 

B.  The  Sequential  Likelihood  Ratio  Test 

The  SPRT  cannot  be  applied  directly  to  test  the  above,  due  to 
the  composite  nature  of  the  hypothesis.  We  will  adopt  as  our  termina¬ 
tion  rule  a  generalization  of  the  SPRT,  the  sequential  likelihood  ratio 
test,  or  SLRT. 


1 .  Two-decision  problem 

We  will  first  consider  the  two-decision  problem 


:  A  =  A  *,  vs. 

H2  :  A  =  -  A  * . 

Let  (X- ..  ,  X  ,  ...,  X  )  be  the  responses  of  n  subjects  administered 
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treatment  1,  and  let  ^22*  ^2n  )  t^ie  responses  of  n^ 

subjects  administered  treatment  2.  We  assume  the  responses  are 
completely  independent,  and  thus  the  joint  likelihood  function  of 


these  variates,  given  n^,  n^  and  0 ,  is  defined  to  be 


n. 


n. 


L  (X;  0)  =  it  f,  (X,4  |  0)  jj  f,  (X 


1  '  li 
i  =  1 


2j 


0) 


j  =  1 


The  likelihood  ratio  for  testing  vs  is  then  given  as 


A  = 


Sup 

A  e  0  e  T 


L  (X;  0  ) 


Sup 

A  E  H2,  0  E  T 


L  (X;  0  ) 


Substituting  (2.1)  and  (2.2)  in  this  gives 

n 


n„ 


A  = 


Sup  “1  “2 

A  e  H1,  0  e  T  TT  f  (X1±  |  0  +  A  /  2)  TT  f  (X2  I  0  -  A  /  2) 

L^J: J  =.  .1 


Sup 


n 


_JL 
i  =  1 


n 


A  e  H  0  e  T  I  I  f  (X.  .  0  +  A  /  2)  I  I  f  (X_  .  |  0  -  A  /  2) 

2  1  1  li  11  2j 


_2 
j  =  1 


Sup  ^1  Tz 

0  e  T  II  f  (X,  .  I  0  +  A*  /  2)  |  |  f  (X0.  |  0  -  A*  /  2) 

II  li  1  II  2 j  1 


A  = 


i  =  1 


j  ..=  1 


Sup  ^1  ^2 

0  e  T  ||  f  (X,  .  |  0  -  A*  /  2)  |  |  f  (X„.  |  0  +  A*  /  2) 

I  I  li  l  l  2j 


i  =  1 


j  =  1 


(2.5) 


Constants  A  and  B  are  chosen  such  that  0<A<1<B<°°.  Sampling 
continues  until  the  first  time  A  i  (A,B),  when  is  accepted  if 
A  <  A  and  is  accepted  if  A  >  B. 


T  1 
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2 .  Three-Decision  Problem 

With  the  addition  of  H  :  A  =  0,  the  test  statistic  becomes 

o 

I 

(A  1»  A  ),  where 

Sup  L  (X;  0  ) 

A  e  H. ,  0  e  T 

A±  =  - - -  ,  i  =  1,2 

1  Sup  L  (X;  0) 

A  e  H  ,  0  e  T 
o 

The  denominator  of  this  expression  is 

Sup  ^1  n2 

0  e  T  |  |  f  (X1:L  |  0)  |  |  f  (X  |  0  ). 

i  =  1  j  =  1 

The  numerator  ofA^  is  identical  to  the  numerator  of  (2.5),  while  the 

numerator  of  A  ^  is  the  denominator  of  (2.5) 

Constants  A  and  B  are  chosen  as  in  the  two-decision  case, 

sampling  continuing  until  max  (a  ,  A0)  t  (A,B).  H  is  accepted  if 

J.  Z  O 

max  (a1,  a9)  <  A*  If  max  (Ax»  A2^  >  B»  ^  is  accepted  if  A  >  and 
is  accepted  if  A^  >  A^. 

A  complete  discussion  of  this  test  may  be  found  in  Chapter  5 
of  Ghosh  (12) . 

It  should  be  noted  that  the  hypotheses  and  could  be 
generalized  to 

H^*  :  A  >_  A*,  vs. 
h2*  :  A  <_  -A  *. 

However,  even  in  the  simplest  cases,  finding  the  suprema  of  the  like¬ 
lihood  functions  under  these  hypotheses  is  very  difficult.  Furthermore, 
termination  rules  become  data-dependent .  The  tests  discussed  in  this 
thesis,  although  based  upon  the  simple  hypotheses,  may  be  considered 
approximations  to  tests  of  the  composite  hypotheses. 


■ 
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C .  Allocation  Rules 

In  this  section,  we  present  the  formulation  of  a  general  data- 

/ 

dependent  allocation  rule  in  discrete  time,  as  given  by  Louis  (15). 


Let  X^,  ...  and  X^*  ***  two  se<luences  °f 

independent  random  variables,  identically  distributed  within  the  sets, 
but  not  necessarily  between  the  sets.  Assume  at  time  n  there  are  n 


X,  ,’s  in  the  sample  and  n_  X_.'s,  so  that  n  =  n,  +  n^, 
li  y  2  2j  12 

of  random  variables  is  defined  as  follows: 

Let 


A  new  sequence 


"n+1  = 


X  if  k  =  1 

1  ,n^+l  n+1 


where 


\x?  >  1 

\  2>n2+i 


1£kn+l  =  0  ' 


ki  =< 


|1,  with  probability  p^  (Z^,  ...»  Z  k  ,  . 


k  ^ 

’  ^i-l 


i,  with  probability  1  -  p.  (Z, ,  ...,  Z.  ,  ,  k  .  ...k 

l  1  l—l  i 


-1) 


i  >  1 


and  p1  =  P, 

0  <  p.  <1,  i  >  1. 

—  i  — 

The  k^'s  constitute  a  sequence  of  indicator  random  variables, 
not  necessarily  either  independent  or  identically  distributed,  which 
specifies  the  order  in  which  treatments  are  assigned.  The  allocation 
probabilities,  the  p^'s,  may  be  allowed  to  depend  upon  both  past  observa¬ 
tions  (Z^’s)  and  past  allocations  (k^’s).  The  allocation  rule  is 
completely  determined  by  the  constant  p  (usually  taken  to  be  1/2),  and 
the  rules  for  calculating  the  p^'s. 

It  is  easily  seen  that  this  formulation  subsumes  both  strict 
alternate  sampling  and  completely  randomized  sampling.  Randomization 
results  from  setting  p^  equal  to  a  constant  for  all  i,  while  strict 
alternation  is  achieved  by  taking 


■ 


. 
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,1  =  1 

l 

The  general  allocation  rule  permits  what  might  be  termed 

"data-dependent  randomization”.  For  administrative  reasons,  this 

possibility  should  be  excluded,  and  we  therefore  require  that 

p.  =  0  or  1  for  i  >  1. 
i 

As  a  general  consideration,  the  p.'s  are  required  to  depend 
upon  the  observations  through  the  likelihood  ratio,  because,  when  a 
test  is  close  to  termination,  we  would  like  to  sample  principally 
according  to  the  "nearest”  hypothesis.  In  other  words,  if  there  is 
strong  evidence  for  either  or  H^,  the  suspected  poorer  treatment 
should  be  assigned  less  frequently  than  the  other.  If  the  trend  is 
towards  acceptance  of  ,  approximate  alternate  sampling  is  desired. 

D.  Performance  Criteria  for  a  Sampling  Plan 

Three  surfaces  which  together  describe  the  relative  merits  of 
a  protocol  may  be  defined,  in  order  to  compose  allocation  rules. 

To  emphasize  the  time-dependence  of  the  likelihood  ratio,  we  shall 
write  A  as  A  (n) .  Define 

N  =  min  {  n  :  A  (n)  t  (A,B)} 

N_^  =  min  {  n  :  A  (n)  i  (A,B)}  ,  i  =  1,2 

These  are,  respectively,  the  total  sample  size  and  the  sample 
sizes  on  each  treatment  at  termination.  Obviously  N  =  +  N^.  The 

average  sample  number  (ASN)  is  then  defined  by 


ASN  =  E  (N). 


■ 
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The  inferior  treatment  number  (ITN)  is  defined  by 

f E  (N.),  where  i  =  <1,  A  <  0 
ITN  =  /  1  L2,  A  >0 

1/2  ASN,  =  0 

Although  we  are  only  considering  two-  and  three-decision 
problems,  we  observe  here  that  the  ASN  and  ITN  are  defined  independently 
of  the  cardinality  of  the  decision  space.  The  third  surface  to  be 
defined,  the  operating  characteristic  (0C),  does  not  have  this  in¬ 
dependence.  We  define  it  here  for  a  k-decision  problem,  with 

hypotheses  H  ,  ...,  H  . 

1  k  k 

Let  S  =  {  (p_ ,  ...,  p.  )  :  Z  p  =  1 ,  p .  >  0,  i  =  1,  ...,  k} 

1  k  i  =  i  1 

Then  0C  :  D  S,  such  that 

OC.  (a)  =  P^  (A)  =  P  (Accepting  H  |  A  e  D  ,  i  =  1,  .  .  . ,  k) 

For  the  two-decision  problem,  this  reduces  to  the  power 
function  of  the  test  and  its  complement,  which  is 

P  (  A  (N)  >_  B  |  A  ) 

In  the  three-decision  problem,  the  OC  is  determined  by  the  pair 

P  (  A  (N)  >_  B  |  A  ,  i  =  1,2) 

The  general  problem  we  are  addressing  is  the  selection  of 
allocation  rules  which  minimize  the  ITN  among  all  allocation  and 
termination  rules  with  equivalent  OC's. 

The  exact  calculation  of  the  operating  characteristic  of  a 
sequential  test  is  rarely  possible.  In  the  case  of  a  SLRT  with  data- 
dependent  allocation,  very  few  advances  have  been  made.  Louis  (15) 
considered  the  problem  of  testing  for  a  difference  in  drift  between 
two  Brownian  motions.  Although  space  does  not  permit  a  full  development 
of  Louis's  results,  we  state  his  principal  result  as  follows: 
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Theorem:  The  SLRT  with  data-dependent  allocation  for  testing  for  a 

difference  in  drift  between  two  Brownian  motions  has  an 

I 

OC  function  given  by  Wald’s  bounds  provided  allocation 
depends  on  the  observations  at  most  through  the  likelihood 
function. 

By  virtue  of  this  theorem,  if  a  test  can  be  shown  to  be 
approximately  equivalent  to  testing  for  a  difference  in  drift  between 
two  Brownian  motions,  the  OC  of  the  test  can  be  approximated  by  Wald's 
bounds . 

E.  Summary 

An  experimental  protocol  for  two-population  sequential  testing 
is  specified  by  three-decision  rules,  of  which  the  admission  rule  is 
non-statistical .  In  this  thesis,  we  restrict  our  attention  to  protocols 
using  sequential  likelihood  ratio  tests  for  termination  rules,  and 
consider  the  choice  of  allocation  rules.  Allocation  rules  may  be 
compared  by  means  of  three  functions:  the  average  sample  number, 
inferior  treatment  number,  and  operating  characteristic. 

In  the  next  two  chapters,  we  present  several  examples  of 


experimental  protocols. 


■ 


' 


Chapter  III  Protocols  for  Selecting 


the  Larger  of  Two  Normal  Means 

I 

In  this  chapter  we  develop  the  SLRT  for  selecting  the  greater 
of  two  normal  means.  Attention  is  given  to  allocation  rules  applicable 
to  two-  and  three-decision  problems. 

A.  The  Two-Decision  Problem 

Let  X, ,  X„ ,  ...»  X  ,  ...  and  Y, ,  Y„ ,  ....  Y  ,  ...  be  two 
1  z  n^  1  z  n^ 

sequences  of  independent  normally  distributed  random  variables  with 
E(X^)  =  y^  and  E(Y_.)  =  P2>  aH  having  variance  1.  Robbins  and 
Siegmund  (16)  have  studied  a  termination  and  allocation  rule  for 
sequential  testing  of  the  hypotheses 


H1 : 

P2>  vs* 

H2  1 

y2 

V 

(3.1) 

(These  results  were  first  announced  by  Flehinger,  Louis,  Robbins  and 
Singer  (11)). 

Reparameterizing  the  observations  in  terms  of 


and  0  =  (\i^  +  y2)/2 

changes  the  hypotheses  (3.1)  to 
:  A  >  0,  vs. 

H2  :  A  <  0. 

1 .  Termination  Rule 

We  present  here  an  alternative  derivation  of  Robbins’  and  Siegmund ’s 
test,  in  the  general  framework  of  Chapter  II.  A  constant  A*  >_  0  is  chosen, 


16 


•  '  -  ' 
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and  the  hypotheses  to  be  tested  become 


:  A  =  A  ,  vs. 

H2  :  A  -  -  A* 


(3.2) 


We  now  develop  the  SLRT  for  testing  vs.  H  .  To  derive  the 
likelihood  ratio,  we  first  find  the  maximum  likelihood  estimator  of  0 
for  A  fixed.  Let  C  =  (2tt)  n^.  Then  the  likelihood  function  is 

nl  n2 

L(X,Y;6,A)  =  C  n  exp[-!s(X  -u  )2]  n  exp[-!j(Y.-u  )2] 

i=l  j=l  J 

nl  n2 

=  C  n  exp[-%(X.-8-A/2)2]  n  exp[-*5(Y.-0+A/2)2] 

i=l  1  j=l  2 


Let  i  =  £n  L 


nl  n2 

-  Hn  C  -  j  I  (X.-6-A/2)2-  i  l  (Y.-0+A/2)2 
i=l  1  j=l  J 


Differentiating  with  respect  to  0,  and  setting  the  result  equal  to 
zero  gives 


n 


n. 


1  2 
l  (X.  -  0  -  A/2)  +  l  (Y.  -  0  +  A/2)  =  0. 


(3.3) 


i=l 


j=l 


,  Ul  x  n2 

Let  X  =  —  T  X.,  and  Y  =  —  Y  Y..  Then  (3.3)  it 

n  u  i  n  u  i 

1  i=l  2  j  =  1  J 

n^X  -  n^0  -  n^  A/2  +  n^Y  -  n^0  +  n ^  A/2  =  0, 


whose  solution  for  0  is 


n  X  +  n0Y  nn  -  n_ 

8  =  1  „  2  -  f  (  1  2) 

n  z  n 


, 
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n  X  +  n  Y  *  n  -  n 

Denote  by  9=  — - - - ■=•  (— - 

+  n  2  n 


«lx  +  V  a*  ,»!  -  "2 
e_  =  - - -  +  -  ( - - - 


\ 


) 


j 


(3.4) 


The  likelihood  ratio  then  becomes 


A(n)  = 


L(X,Y;0+,A  ) 
L(X,Y;0  ,-*) 


ni  *  * 

=  n  exP{-Js[(x.-e  -  f  )2  -  (x.-e_+  £  )2]} 

i=l  z  1  -  Z 


n2  ■k  k 

n  exp{-%[  (Y  -§++  |  )2  -  (Y.-0_-  |  )2]} 
j=l  J  J 


(3.5) 


Upon  substitution  for  0+  and  0_,  the  first  product  in  (3.5)  becomes 

nl  n  X  +  n  Y  n  -  n  * 

n  exp{  (— - —  -  X.  )  (A  (— - -)  -  A  )} 

i-1  n  i 

T_  jJ.  _ 

1  -2A  n  n  X  +  n  Y 

=  n  exp{ - - —  ( - - - X  )} 

i=l  n  n  i 


-2A  n  n  n  X  +  n0Y  n_n0 

r  ±  Z  /  L  ^  \  I  n  »  A  -1-  ^  Vi 

=  exp  [ -  ( - )  +  2A  -  X] 


n 


n 


n 


(3.6) 


Similarly,  the  second  product  in  (3.5)  becomes 


2A  n1 n  n.,X  +  n0Y  .  nnn„ 

r  ±  Z  /  L  ^  \  r|A  J-  A  y 

exp  [ - - -  ( - - - )  -  2A  — —  Y] 


n 


n 


n 


(3.7) 


Multiplying  (3.6)  and  (3.7)  gives  the  likelihood  ratio 


k  nin2  —  — 

A(n)  =  exp[2A  — —  (X  -  Y)]. 


n 


The  continuation  region  of  the  test  is  given  as  A  e  (A,B),  where 
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0  <  A  <  1  <  B  <  °°.  If  we  define 

£n  A  ,  ,  £n  B 
a  =  - —  and  b  =  - —  , 

2A  2A 

the  test  may  be  based  upon  the  statistic 

nin2  —  — 

C(n)  =  — —  (X  -  Y)  ,  (3.8) 

with  termination  occurring  the  first  time  £(n)  £  (a,b) .  If  £(n)  <_  a, 
is  accepted.  If  £ (n)  >_  b,  is  accepted. 


2 .  The  Operating  Characteristic  of  the  Test 

The  statistics  £(n)  are  linear  functions  of  jointly  normal 
random  variables,  and  are  therefore  jointly  normally  distributed  for  all 
n.  We  have 


nln2  -  - 

E  r(n)  =  E(X  -  Y) 

A  n 


nln2 

n 


A  ,  and 


(3.8) 


VarAa(n))  =  -±f  -  Var(X  -  Y) 

n 

2  2 

nln2  —  — 

=  — ^  [Var (X)  +  Var(Y)] 

n 


2  2 

nln2 

2 

n 


nln2 


n 

Thus,  the  distribution  of  c(n)  does  not  depend  on  0,  and  the  test  is,  in 
fact,  a  SPRT  for  testing  vs.  H2 .  (This  latter  fact  is  perhaps  more 


v 


. 
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clearly  seen  via  the  derivation  by  Robbins  and  Siegmund  (16) ,  who  employ 
a  linear  transformation  of  the  observations  to  eliminate  0.) 

I 

Therefore,  for  any  allocation  rule  for  which  termination  occurs 
with  probability  1,  Wald’s  bounds  approximate  the  error  probabilities  of 
the  test.  Thus,  we  have 

P_A  (accept  H^)  <  B  ^ [ 1  -  (accept  H^)] 
and  P^  (accept  H^)  <  A[1  -  P_A  (accept  H^)] . 


There  is  approximate  equality  in  these  statements,  so  simultaneous 
solution  gives 


P_A  (accept  H1) 


and 


(accept  H^) 


^  A(B  -  1) 
B  -  A 


In  terms  of  a  and  b,  these  become 


_  ,  TT\^  1  -  e 

P_A  (accept  H1)  =  - * - 


* 

2A  a 


*  “IT  ’  and 

2A  b  2A  a 


-  e 


2A*a  2A*b  _  . 

PA  (accept  H2)  =  6  6 


* 

2A  b  2A  a 


(3.9) 


(3.10) 


e  -  e 

These  errors  probabilities  do  not  depend  on  the  particular 
allocation  rule  used,  provided  the  rule  depends  on  the  observations 
at  most  through  the  differences  X^  -  Y  . 

If  the  allocation  rule  is  symmetric  in  and  Y  ,  and  a  =-b,  we 
have  P_A  (accept  H^)  =  PA  (accept  F^) .  Then  (3.9)  and  (3.10)  give 

P  (error)  <  - ,  A  ^  0.  (3.11) 

A  i  *  -2Aa 

1  +  e 


. 
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Again,  there  is  approximate  equality. 

3 .  Estimates  of  Expected  Sample  Sizes 
If  F 


denotes  the  a-algebra  generated  by  the  x-r’s 
n  ,n_  1 

1  2  n,n„ 


and  the 


y^’s,  it  may  be  shown  that  the  observed  process  (c(n)  -  n 


12 


A,  F 


nl’n2 


is  a  martingale.  Robbins  and  Siegmund  exploit  this  fact  to  estimate 
expected  sample  sizes  in  the  symmetric  case  referred  to  in  the  last 
paragraph  of  the  preceding  section. 

Their  results  can  be  summarized  as  follows: 

1.  For  A  i  0 


N..1L  ,  2bA 

F  (  1  2  s  ^  b  e  -  1. 

AV  N  }  A  ^  2bA 

e  +1 


(3.12) 


If  A  =  0 


V-^r’ 1  b' 

2.  Strict  alternate  sampling  has  N  =  N0  =  N/2,  so  (3.12)  gives 


,,  2bA  . 

VN)  -r  (^2biTf) 

e  +1 


2bA  -  1 

and  EA(Nl)=  Ea(N2)  -  - 

e  +1 


(3.13) 


3.  Since  min(N. ,N0)  > 


N1N2 


1*  V  -  N 


,  for  any  allocation  rule  we  have 


N  N 

min[Ei(N1),EA(N2)]  >  E&  (-5-) 


N 


(3.14) 


4.  Simple  calculations  show  that,  for  all  -  N 


. 


Therefore 


(3.15) 


Equality  holds  only  if  P  (N.^^)  =  ^ * 

There  are  two  important  consequences  of  these  results.  First, 
(3.12)  through  (3.14)  show  that  the  expected  sample  size  on  each  treatment 
must  exceed  approximately  one-half  the  expected  number  on  that  treatment 
under  alternate  sampling.  Second,  (3.15)  states  that  the  ASN  is 
(approximately)  minimized  under  alternate  sampling. 

4 .  Allocation  Rule 

Robbins  and  Siegmund  studied  an  allocation  for  the  symmetric  test 
discussed  in  the  preceding  sections.  We  formulate  this  rule  in  the 
general  form  of  Chapter  II.  A  constant  c  >  0  is  chosen. 


Define  I(S) 


Let 


probability  I[ 


where 


0,  with  probability  1  —  I [ 


1,  i  =  1 


\  0,  i  =  2 


That  is,  we  begin  by  observing  and  Y^.  For  succeeding  observations 
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if  termination  has  not  occurred  we  observe  an  X  if 

n,  -  n_  *  -v 

_1 2  C(n) 

n  —  c  * 

otherwise,  a  Y  is  observed. 

Monte  Carlo  results  on  sample  sizes  and  OC's  were  obtained  by 
Robbins  and  Siegmund  (16) .  A  selection  of  these  results  is  presented  in 
Table  1  (see  Table  1,  page  24).  The  sampling  rule  appears  to  come  very 
close  to  achieving  the  theoretical  bounds  given  in  section  A. 4. 

B .  The  Three-Decision  Problem 

We  now  add  the  "indif f erence"  hypothesis  H  :  A  =  0  to  the 

o 

problem  stated  in  section  A.l.  Flehinger  and  Louis  (9)  have  studied 
a  termination  rule  and  a  class  of  allocation  rules  for  the  three- 
decision  case. 


1.  Termination  Rule 


We  shall  use  the  same  nomenclature  as  in  the  two-decision 


problem.  Denoting  by 


0  = 
o 


n  X  +  n2Y 
n 


the  likelihood  ratios  are  this 


A-^(n)  « 


L(X,Y;e+,A  ) 
L(X,Y;0o,O) 


A2 (n)  = 


L(X,Y;§  ,-A  ) 


L(X,Y;6o,0) 


9 


and 


, 
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TABLE  1 

MONTE  CARLO  RESULTS  FOR  ROBB INS -SIEGMUND 
ALLOCATION  RULE:  SYMMETRIC  TEST  WITH  b  =  6 

A  OC  ASN  ITN 


WA(1) 

c=6 

c=7 . 2 

WA 

(alternate) 

o 

n 

ON 

c=7 . 2 

WA 

(alternate) 

c=6 

CM 

• 

r^. 

II 

u 

0 

- 

- 

- 

144 

216 

191 

72 

108 

95 

.1 

.77 

00 

1-"- 

• 

.75 

128 

197 

166 

64 

86 

80 

.25 

.95 

.96 

.97 

86 

1 

146 

126 

43 

34 

36 

.5 

.50 

.996 

1.00 

48 

109 

78 

24 

16 

16 

1.0 

1.00 

1.00 

1.00 

!  24 

| 

71.1 

45.2 

12 

7.2 

7.7 

(1)  WA  =  Wald  approximation 
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A1(n)  = 


L(X,Y;6+,A  ) 
L(X,Y;6o,0) 


L(X,Y;0+,A*) 


n 


*  n2  * 

C  n  exp[-Js(X.-0  -  J  )2]  n  exp[-Jg(Y.-e  +  |  )] 

•  f  1  •  ^  •  1  "1  '  ^ 

1=1  J=1 


9+  +  2  = 


*  n,X  +  n„Y  A*  n,  -  n?  * 

4-  — 

2 


-  1) 


9+  -  2  = 


n. 


1 

z 

a  .  i  z 

n 

2  (  n  } 

n  X 

+  n„Y 

*  n„  -  n„ 

2 

A  ,  1  2 

n 

2  (  n 

_ 

* 

n  X 

+  n2Y 

+  n2A 

n 

JL 

A 

nX 

+  n_Y 

-  n.  A 

1 

2 

1 

ni 

n 

n  X  +  n  Y 

(X1  '  (  n 

n  X  +  n?Y  , 
(  1  Z 

n  i 

exp [-^ 

=i 

ovn  f  (  T  V 

n. 


i=l 


n 


n  X  +  n  Y 

))"]  n  expHsCY.  -  (— - —))Z] 

j=l  J 


n  X  +  n  Y  +  n9A  9 

-  (— - f - — )]  )} 

l  n 


n 


n  X  +  n  Y  n  X  +  n?Y  -  n  A  9 

n  {3s([Y,  -  (— - — )T  -  [Y.  -  (-± - - - — )T} 

J  n  j  n 


j=l 


The  first  product  equals: 


n. 


n  exp{^(2X.  -  2(— - - 

i-1 


n,X  +  n„Y  n  A  n  A 

-)  -  ~— )  (^— )  > 
n  n 


(3.16) 


The  second  product  equals: 


n„ 


II  exp{^(2Y  -  2  ( - — 

j=l  3 


jj, 

n, X  +  n„Y  n^A  n  A 


-)  +  -=—)  (~  -z—) } 


n  n 


(3.17) 
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(3.16) 


(3.17) 


Thus , 


Because 


chosen. 


becomes 


exp{ 


*  n 

n2A  1 


n 


n,X  +  n^  n2A 


l  [X,  -  "  )  -  -A— ]} 


i=l 


n 


2n 


=  exp{ 


n2A  -  nl  -  -  nln?A 

—  [nlX  "  V  (nlX  +  n2Y>  -  ^25“  1} 


"2A  rnln2-  nln2-  nlV  ,, 

exp{ -  [ -  X - Y - - -  ]  } 


n  n 


n 


2n 


nln2  *  —  —  a* 

exp{  2~  A  [X  -  Y  -  -  ]  } 


n 


becomes 


*  n 
-n,A  2 


exp{- 


1 A  2  n  _  _  n  A 

—  l.l  [Yj  --  (nix  +  n2Y)  +-ir]} 

3=1 


2  * 

-ll  n  A  _ 

exp{ - ;; —  [Y 


n 


X  +  |  ]} 


n  n  ^  _  _  * 

A  (n)  =  exp{ - -  A  (X  -  Y  -  —  )} 

In  z 

A1(n) 

- — 7— r-  =  A(n)  (from  the  definition),  we  have  immediately 
A2(n) 

n  n  ^  * 

A2(n)  =  exp{— —  A  (Y  -  X  -  —  )} 

As  in  the  two-decision  problem,  two  constants,  A  and  B,  are 


such  that  0<A<  1<  B<°°.  The  procedure  terminates  when 


max(A^,A2)  leaves  the  interval  (A,B). 


' 


27 


2.  A  Class  of  Allocation  Rules 


We  now  define  the  rules  R  .  A  constant  y  is  selected  with 

Y 


0  £  y  £ 

1. 

- 

Define  A  =  X  -  Y 

ni 

• 

n 

2 

Let 

-  < 

1, 

n1  -  n2  _>  y  (n  +  1) 

or 

ni 

Pn+1 

0, 

n1  -  n2  £  -  y(n  +  1) 

or 

ni 

where  p  , , 
n+1 

=  Pn+l(Zl . Zn’kl’ 

• • • ,kn) 

2 

l2 


When  y  =  0,  the  rule  is  strict  alternation,  and  when  y  =  1,  it 
assigns  the  leading  treatment  at  every  arrival. 

Flehinger  and  Louis  (9)  obtained  Monte  Carlo  results  for  the  R 
rules,  using  the  stopping  rule  discussed  in  the  preceding  section. 

Table  2  (see  Table  2  on  page  28)  contains  a  selection  of  their  results. 


Y 


. 
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TABLE  2 

MONTE  CARLO  RESULTS 

NORMAL  DISTRIBUTION  -  THREE-DECISION  PROBLEM 
WITH  A  =  0.1,  B  =  30 


A  OC  ASN  ITN 


y 

.2 

.5 

.2 

.5 

.2 

.5 

0 

.05 

.05 

127 

160 

64 

80 

.125 

.13 

.14 

141 

180 

63 

66 

.25 

.43 

.43 

164 

211 

68 

62 

.50 

.94 

.94 

107 

136 

43 

36 

1.0 

1.00 

1.00 

40 

51 

16 

13 

CHAPTER  IV  PROTOCOLS  FOR  SELECTING  THE  BETTER  OF 


TWO  EXPONENTIAL  PARAMETERS 

/ 


In  this  chapter  we  examine  allocation  and  termination  rules  for 
selecting  the  better  of  two  treatments,  when  patients’  survival  times  are 
exponentially  distributed.  Our  attention  shall  be  given  entirely  to  the 
three-decision  problem,  and  the  comparison  of  several  allocation  rules. 


A.  The  Testing  Situation 

We  consider  a  situation  in  which  suitable  patients  arrive  for 
treatment  at  a  rate  of  one  patient  per  unit  time.  There  is  no  loss  of 
generality  in  this  assumption,  since  time  units  may  be  re-defined  as 
necessary.  Patients  are  assigned  treatments  by  an  allocation  rule,  and 
their  survival  times  after  treatment  are  recorded. 

We  assume  survival  time  of  patients  receiving  treatment  i 
follows  an  exponential  distribution  with  instantaneous  death  rate  A.. 
That  is,  if  X  is  survival  time  then  the  density  of  X  is 


f(x)  = 


-A  .x 

A .  e  1  ,  x  >  0 
1  — 


0,  x  <  0 


We  will  also  have  occasion  to  use  the  parameter  vu 
mean  of  the  distribution. 


A.\  which  is  the 

l 


As  in  previous  chapters,  the  subjects  are  assumed  to  be 


independent . 


A  constant  p  >  1  is  chosen  which  indicates  a  medically 


significant  difference  in  lifetime.  We  wish  to  test  the  three  hypotheses 


29 


■ 


30 


(4.1) 


(These  hypotheses  can  be  written  in  the  form  of  section  2. A,  but 
the  formulation  above  is  more  natural  for  the  exponential  case.) 

It  is  further  assumed  that  data  on  each  patient  is  updated 
continuously.  This  being  the  case,  at  each  time  n,  we  have  the  following 
data: 


n  ,  i  =  1,2,  the  number  of  subjects  assigned  treatment  i, 

with  nn  +  n_  =  n; 

In  2n 

d^^,  i  =  1,2,  the  number  of  deceased  subjects  on  treatment  i; 


X..  ,  i  =  1,2,;  j  =  1,2, ...,n.  ,  the  survival  time  of  the  j 
ijn  m 


th 


patient  on  treatment  i  up  to  time  n 


(4.2) 


To  distinguish  between  active  patients  and  those  who  have  died, 


the  X..  ’s  are  separated  into  two  subsets  for  computational  purposes: 
ijn 


x. .  ,  i  =  1,2;  j  =  1, . . . ,d.  ,  the  survival  time  after 

ijn  J  in’ 

treatment  of  a  deceased  patient; 

y.  .  ,  i  =  1,2;  j  =  l,...,n.  -  d.  ,  the  time  since  treatment 

Jijn*  *  m  in’ 

of  an  active  patient  (4.3) 

The  survival  times  are  summarized  by  the  statistics  t^n  and  t^^^, 

where 


t. 

m 


n . 
in 


I  X. .  ,  i 
j=l  ljn 


1,2 


(4.4) 


■ 
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That  is,  t.  is  the  total  survival  time  on  treatment  i  up  to 
m  r 

time  n.  In  the  sequel,  we  shall  suppress  the  n  subscripts  for  simplicity 

I 

of  notations. 

1 .  Graphical  Representation  of  Data 

There  is  a  convenient  way  of  representing  this  type  of  data 
graphically.  Figure  1  illustrates  this  method  for  a  small  set  of  data, 

12  observations  with  y  =  10,  and  11  with  y  =  5,  under  alternating 
treatment  allocation  (see  Figure  1  on  page  32) . 

The  sample  from  y  =  10  gave  values  4.3,  4.7,  4.7,  etc.,  and  the 
sample  from  y  =  5  gave  values  4.5,  5.0,  1.4,  etc.  At  time  24,  shown  on 
the  graph,  there  are  9  patients  surviving  from  the  original  23;  4  of  the 
survivors  had  received  treatment  2,  the  rest  treatment  1. 

B.  Maximum  Likelihood  Estimation  of  the  Death  Rate 

A  maximum  likelihood  estimate  of  the  instantaneous  death  rate 
is  required  in  the  derivation  of  the  SLRT.  The  "usual"  estimate  is 
not  applicable  in  this  situation,  since  the  data  are  partially  censored 
by  survival. 

This  estimation  problem  was  studied  by  Bartholomew  (5) ,  whose 
main  result  is  expressed  in  the  following 

Theorem:  Let  X^. ,  j  =  l,...,n  be  a  random  sample  from  an  exponential 
distribution  with  death  rate  X.  Define  d  and  t  as  in  (4.2) 
and  (4.4).  Then  the  maximum  likelihood  estimate  of  X  is 

X  =  d/t,  (4.5) 

provided  entry  times  are  independent  of  previous  data. 


< 
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/ 

Graphical  Representation  of  Survivor  Data 


15.0 

13.0 

12.0 


6.0 

5.0 

4.0 

3.0 

2.0 

1.0 


Time 


- ►Surviving  Subject 

jLt1  (odd  numbers)  =  10 

jn2  (even  numbers)  =  5 

n1  =  12  d-,  =  7  t1  =  59.0 

n2  =  1 1  d2  —  7  t2  =  49.7 
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Proof :  Define  P  -  P [subject  has  died  |  entry  X.  ago],  j  =  l,...,n 

J  3 

-AX. 


Q.  =  1  -  P.  =  e  J 

J  3 


^  ^  * )  j  1  j  •  •  •  j  n 

J  J 


e .  = 
J 


t°’  X3  =y3 

where  x.  and  y.  are  defined  as  in  4.3. 


n 

Note  that  Y  e .  =  d. 

j-i  J 

Under  the  assumption  of  independence,  the  likelihood  function  may  be 
written: 


n  e .  1-e .  -AX .  e . 
L(A)  =  n  P.JQ.  J{Ae  J/P.}  3 
2=1  3  3  3 

n  -Ay.  1-c .  -Ax.  e. 


=  n  (e  J)  3  (Ae  J) 

j=l 


Jn  J 


n 


n 


=  A  exp [-A (  l  (1  -  e.)y.  +  £  e.x.)] 


j=l  3  3  j=l  3  3 


L(A)  =  Ade  At 


Let  £(A)  =  £n  L(A)  =  d  £n  A  -  At 


Then  If  =  i  -  E 


Setting  this  equal  to  zero  gives  the  desired  result. 


(4.6) 


C .  The  SLRT  for  the  Exponential  Parameters 

Returning  to  the  setting  of  sections  4. A,  we  now  develop  the 
SLRT  for  the  hypotheses  (4.1).  Reinstating  the  i  subscripts  in  (4.6) 
gives 

d.  -A . t . 

W  =  h*6  1  1  »  1  =  i*2 
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Assuming  that  treatment  assignment  is  independent  of  past  data,  the  joint 
likelihood  function  of  A^  and  A^  is 

/ 

L(X1,X2)  =  l1(x1)l2(x2) 

dl  d2 

=X^  exp(-X1t1)X  exp(-X  t  ) 

When  A^  =  pA^,  this  becomes 
^2  ^1+C^2 

L  =  p  \1  exp[-A1(t1  +  pt2)] 

For  fixed  d^,  t^,  this  is  maximized  by 


A,  = 


dl  +  d2 


1  tl  +  pt2 


Recalling  that  we  are  interested  in  the  hypotheses  (4.1),  we  must 

■k  k 

find  the  maximum  values  of  L(A^,A2)  subject  to  p  =  1,  p  =  p  ,  p  =  1/p  , 
respectively.  These  maximum  values  are 


d  +  d  d  +d 

Lo  =  (T~T^  exp("di  '  d2) 


d0  d  +  d  d  +d 

L.  =  p*  ( - * — )  expC-dj^  -  d2) 

h  +  p  l2 


d  d  +  d0  d  +d 

L2  =  p*  1(~z - )  exp(-d1  -  d2) 

P  si  +  t2 

The  likelihood  ratio  pair  (A  ,A2>  is  thus  given  by 


L_  d0  t-  +  t„  dn+dn 

.  1  *2,1  2  s  1  Z 

Ai  =  TT  =  p  ( - - } 

1  o  t,  +  p  t, 


' 
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d  t  +  t  d  +d 
*1/1  l  N  1  L 
P  (r~z - ) 

P\  +  t2 


l 


As  in  section  2.B.  we  terminate  sampling  as  soon  as  max(K^yA^) 
leaves  the  interval  (A,B)  where  0<A<  1<  B<°°.  This  test  is  due  to 
Flehinger  and  Louis  (10),  and  is  used  in  conjunction  with  each  allocation 
rule  in  the  next  section. 


D.  Allocation  Rules 


1 .  The  Ry  Rules 

A  class  of  allocation  rules,  proposed  and  examined  by  Flehinger 

and  Louis  (10) ,  are  the  analogues  of  the  R^  rules  for  the  normal  case 

described  in  the  previous  chapter.  As  in  that  case,  a  constant 

0  <_  y  <_  1  is  chosen  and  kept  fixed  throughout  the  experiment.  At  the 

st 

time  of  arrival  of  the  (n  +  1)  subject  a  treatment  is  assigned  to 
him  according  to  the  probability  Pn+1  of  assigning  treatment  1: 

1,  d2  ~  di  —  +  1)  »  or  I  dx  “  d2 1  <  "Y (n  +  1)  and  £  ^2 

0,  dx  ”  d2  £  y(n  +  1) >  or  ldx  ”  d2 I  <  Y^n  +  and  ^2  —  ^1 

1 

P1  =  2 

The  effect  of  this  rule  is  to  assign  the  treatment  on  which  fewer 
subjects  have  died,  if  the  difference  between  death  numbers  is  sufficiently 
large.  If  the  difference  is  small,  the  rule  chooses  the  treatment  with 
the  lower  estimated  death  rate. 

The  indeterminacy  in  the  case  d1  =  d2  =  0  is  avoided  by  assigning 
alternately  until  at  least  one  death  has  been  observed.  Equivalently, 


■ 


. 
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~  -1 

we  may  base  the  rule  on  =  A^  ,  taking  Bartholomew's  (5)  version 
of  this  estimate  in  the  case  d.  =0.  That  is,  if  d.  =  0,  then 

l  l 

VU  =  t  ,  i  *  1,2. 

a.  Characteristics  of  the  R  rules 

Y 

Estimated  operating  characteristic  (OC) ,  average  sample 
number  (ASN) ,  and  inferior  treatment  number  (ITN)  surfaces  were  given 
by  Flehinger  and  Louis  (10),  results  which  we  were  able  to  confirm  with  a 
reasonable  degree  of  accuracy  with  our  own  Monte  Carlo  trials.  Detailed 
output  from  these  trials  revealed  several  noteworthy  features. 

(i)  The  ITN  appears  to  be  minimized  for  some  y  near  0.5,  for  all 
A^  and  X^. 

(ii)  For  all  A ^  and  X^,  the  ASN  increases  with  y.  Analogous  to  the 
normal  theory  results  of  Robbins  and  Siegmund,  the  ASN  seems  to 
be  minimized  by  strict  alternate  sampling.  (See  Table  3  on  page  37.) 

(iii)  Sampling  tends  to  proceed  in  long  runs  of  each  treatment. 

(iv)  Let  p  =  ^^^2  ^xec**  Then  the  ratio  ITN/ASN  decreases  in  X ^ 
for  a  fixed  p. 

(v)  For  a  given  stopping  rule,  the  OC  function  depends  only  on  p. 

(See  Table  4  on  page  38.) 

In  summary,  it  appears  that  the  R^  rules  oversample  on  the  better 
treatment,  particularly  for  small  values  of  X and  A^.  This  may  be 
due  to  the  tendency  to  stick  on  a  single  treatment.  If  the  sampling 
pattern  were  more  mixed,  the  amount  of  information  on  each  treatment 
would  increase  at  all  times.  The  choice  of  y  presents  a  considerable 
problem.  Increasing  y  from  0.2  to  0.5  gives  a  small  decrease  in  the 


.1  *  t  *  >  - 
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TABLE  3 

AVERAGE  SAMPLE  NUMBERS  UNDER  VARIOUS 
ALLOCATION  RULES,  FOR  SELECTED  MEAN  LIFETIMES 


h1 

A'1 

Rule 

ASN 

10 

10 

Alternation 

75. 

Ro 

77. 

R.  2 

78. 

R.5 

98. 

10 

20 

Alternation 

68. 

Ro 

69. 

R.  2 

78. 

R.  5 

108. 

100 

200 

Alternation 

136. 

Ro 

137. 

R.  2 

183. 

R.5 

310. 

100 

50 

Alternation 

111. 

Ro 

110. 

R.  2 

127. 

R.  5 

199. 

Termination  by 

SLRT  with 

A  =  0.1,  B  =  30,  p  =  2 

1 

8 

7 

0 

9 

5 

1 

8 

3 

4 

0 

3 

0 

1 

3 

7 


2 
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TABLE  4 

OPERATING  CHARACTERISTICS  OF  SELECTED  ALLOCATION  RULES, 

FOR  SELECTED  MEAN  LIFETIMES 


p 

Rule 

Probability  of 
A~*  =  10 

Selecting  Better  Treatment 

Aj1  =  50 

3 

Alternation 

1.00 

1.00 

R.2 

1.00 

1.00 

R.5 

1.00 

1.00 

2 

Alternation 

0.94 

0.98 

R.  2 

0.96 

0.91 

R.  5 

0.92 

0.91 

1 

Alternation 

(0.95) 

(0.98) 

R.  2 

(0.96) 

(0.96) 

R.  5 

(0.96) 

(0.96) 

1/2 

Alternation 

0.95 

0.95 

R.  2 

0.93 

0.93 

R.5 

0.97 

0.96 

1/3 

Alternation 

1.00 

1.00 

R.  2 

1.00 

1.00 

R.  5 

1.00 

1.00 

1.5 

Alternation 

0.52 

0.47 

R.2 

0.53 

0.62 

R.  5 

0.48 

0.58 

2/3 

Alternation 

0.58 

0.52 

R.2 

0.48 

0.54 

R.  5 

0.55 

0.51 

Termination  by  SLRT  with  A  =  0.1,  B=30,  p  =2 
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ITN  at  the  expense  of  a  (perhaps  unacceptable)  increase  in  the  ASN. 

We  now  examine  two  rules  which  avoid  the  use  of  an  arbitrary 

l 

constant . 

2 .  Likelihood-Bias  Sampling 

Hoel  and  Weiss  (14)  used  likelihood-bias  (L-B)  sampling  in 

& 

conjunction  with  a  modification  of  procedure  P  of  Bechhofer,  Kiefer, 

D 

and  Sobel  (6).  In  the  two-population  case,  this  procedure  calculates 
P^  and  after  each  observation,  where 

P^  =  P[correctly  selecting  treatment  i  |  observations],  i  =  1,2. 
L-B  sampling  uses  the  posterior  probabilities  P^  and  P^  to 
adjust  the  sampling  ratio  by  defining 
P(assigning  treatment  1) 

-  V(pi  +  V- 

Note  that  this  is  a  randomized  allocation  rule. 

Using  the  SLRT,  Sobel  and  Weiss  carried  out  Monte  Carlo  trials  to 
compare  L-B  sampling  with  the  rules.  They  found  the  R^  rules  to  be 
markedly  superior  in  terms  of  both  the  ASN  and  the  ITN. 

3.  Rule  R 

e 

This  rule  attempts  to  "mix"  the  observations  better  than  the  R^ 
rules,  while  avoiding  the  use  of  an  arbitrary  constant.  Unlike  other 
allocation  rules,  it  does  not  depend  on  the  observations  through  the 
likelihood  ratio.  Instead,  it  forces  the  sampling  ratio  n^n  to 
approach  the  ratio  y  /(I^  +  y2) .  This  was  motivated  by  the  observation 
during  the  initial  Monte  Carlo  trials  on  R  ,  that  these  two  ratios  tended 
to  be  approximately  equal. 


. 


40 


Define  a„  = 


a„  = 


nl  +  1 


n  +  1 


n. 


/N  ,  /V 

yl  +  y2 


n  +  1  U1  +  P2 


Then  p  -  1/2  and  p  ...  =  1  -  p  ,  if  either  of  d,  or  d„  is  zero, 
1  n+1  n  12  9 


and 


=  f1’  “l  <  a2 

yn+l  jo.  ai  >  a2 


if  neither  d^  nor  d^  are  zero. 

Monte  Carlo  trials  indicated  that  this  rule  succeeds  well  in 
mixing  the  observations.  However,  when  the  first  subject  dies  almost 
immediately,  y^/ (y  +  y  )  is  close  to  0  or  1,  and  essentially  no  mixing 
occurs  at  all;  the  procedure  may  in  fact  fail  to  terminate. 

The  problem  is  due  to  the  instability  of  y^/ (y^  +  y^)  for  small 
sample  sizes.  We  can  ensure  greater  stability  by  alternating  until  a 
sufficient  number  of  deaths  has  been  recorded;  that  is,  alternation  ceases 
when 


n  =  min{n  :  max(d,  ,d_  )  =  k} . 

In  2n 

k  is  chosen  so  that  the  procedure  cannot  terminate  under  alternation,  by 
solving  A (n)  <  B  for  k.  This  gives 
k  =  [log  .  B]  +  1, 


the  square  bracket  denoting  the  greatest  integer  function. 

Monte  Carlo  results  incorporating  this  "delay"  showed  that  rule 
Rp  gives  a  lower  ASN  than  the  R^  rules,  but  the  R^  rules  tend  to  have 


Y 


lower  ITN ’ s . 


, 
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4.  Delayed  R  Rules 

1 

Introducing  a  delay  to  rule  R  gave  a  substantial  reduction  in 

/ 

ASN.  It  is  therefore  reasonable  to  ask  if  such  a  delay  would  similarly 

improve  the  R  rules.  Monte  Carlo  trials  showed  that  this  was  the  case, 

Y 

along  with  a  reduction  in  the  ITN. 

Tables  5,  6,  and  7  (shown  on  pages  42,  43,  and  44)  give  Monte  Carlo 

results  for  rules  R  ,  R  ,  and  R  9,  all  with  and  without  delays.  In  all  cases 

6  U  •  Z 

: k 

p  =2,  and  termination  took  place  with  A  =  0.1,  B  =  30.  In  delayed  trials, 
the  delay  used  was  5  deaths.  The  delayed  R^  rules  have  lower  ASN’s  and 
ITN’s,  with  relatively  little  change  in  their  operating  characteristics 
from  the  non-delayed  rules.  Rule  R^  is  competitive,  especially  for  values 
of  p  greater  than  2  or  less  than  1/2. 

E.  Summary 

It  is  possible  to  reduce  the  inferior  treatment  number  of  a  two- 

population  sequential  test,  by  using  data-dependent  treatment  allocation. 

In  the  case  of  comparing  the  parameters  of  two  exponential  populations,  we 

have  presented  two  rules,  delayed  R  ,  and  delayed  R  ,  which  improve  upon 

e  y 

previously  published  results. 
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TABLE  5 

CHARACTERISTICS  OF  RULE  R  WITH  AND  WITHOUT  DELAY 


,,  -1  - 
9 

200 

ASN 

4 

V1 

P 

Delayed 

Not 

Delayed 

600 

3 

148 

169 

500 

2.5 

156 

172 

400 

2 

198 

213 

300 

1.5 

223 

231 

266.67 

1.333 

221 

223 

200 

1 

187 

182 

150 

.75 

202 

187 

133.33 

.667 

191 

197 

100 

.5 

138 

134 

80 

.4 

109 

114 

66.667 

.333 

92 

94 

o 

i — i 

ii 

i — i 

i 

CM 

30 

3 

55 

67 

25 

2.5 

60 

76 

20 

2 

80 

90 

15 

1.5 

108 

102 

13.333 

1.333 

106 

106 

10 

1 

80 

87 

7.5 

.75 

90 

95 

6.667 

.667 

97 

96 

5 

.5 

71 

77 

5 

.4 

56 

67 

3.333 

.333 

51 

55 

ITN  OC 


Delayed 

Not 

Delayed 

Delayed 

Not 

Delayed 

44 

48 

1.00 

1.00 

48 

51 

.99 

1.00 

79* 

81 

.96 

.94 

91 

100 

.48 

.54 

96 

97 

.30,  . 

1  .37 

- 

- 

.04 

.06 

88 

81 

.22 

.22 

77 

79 

.58 

.55 

47 

45 

.96 

.92 

34* 

37 

1.00 

1.00 

27 

27 

1.00 

1.00 

15 

18 

1.00 

1.00 

18* 

21 

1.99 

.99 

27* 

31 

.96 

.92 

43 

40 

.62 

.60 

45 

46 

.35 

.36 

- 

- 

.02 

.07 

38 

41 

.29 

.33 

39 

44 

.51 

.56 

24 

24 

.95 

.93 

16* 

18 

1.00 

.99 

13 

13 

1.00 

1.00 

Termination  by  SLRT  with  A  =  0.1,  B  =  30. 
100  trials  per  parameter  pair 

•Jf 

significant  at  p  <  0.05 
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TABLE  6 

CHARACTERISTICS  OF  RULE  R  WITH  AND  WITHOUT  DELAY 

0 


II 

I- 1 

1 

CN 

«-< 

200 

ASN 

ITN 

OC 

A  -1 

1 

P 

Delayed 

Not 

Delayed 

Delayed 

Not 

Delayed 

Delayed 

Not 

Delayed 

600 

3 

148 

155 

43* 

48 

1.000 

.999 

500 

2.5 

157 

171 

51* 

58 

.995 

.993 

400 

2 

194 

196 

73* 

76 

.925 

.949 

300 

1.5 

226 

231 

100 

102 

.525 

.573 

266.67 

1.333 

222 

222 

103 

103 

.345 

.307 

200 

1 

182 

185 

- 

- 

.045 

.053 

150 

.75 

189 

198 

90 

91 

.265, 

.302 

133.33 

.667 

189 

196 

83* 

86 

.005 

.580 

.570 

100 

.5 

142 

146 

55* 

58 

.920 

.950 

80 

.4 

113 

113 

39 

41 

.995 

.994 

66.67 

.333 

92 

97 

29* 

33 

1.000 

.999 

X2 

10 

30 

3 

48 

47 

19 

19 

1.000 

.999 

25 

2.5 

56 

55 

24 

24 

.990 

.993 

20 

2 

68 

71 

31 

32 

.945 

.943 

15 

1.5 

100 

99 

49 

48 

.565 

.572 

13.333 

1.333 

93 

97 

46 

48 

.340, 

.005 

.298 

10 

1 

76 

76 

- 

- 

.070 

.050 

7.5 

.75 

90 

94 

44 

46 

.305 

.247 

6.667 

.667 

96 

96 

47 

47 

.555 

.573 

5 

.5 

65 

65 

30 

31 

.93 

.947 

4 

.4 

45 

50 

20* 

23 

.995 

.995 

3.333 

.333 

36 

45 

16 

21 

.995 

.998 

0.1,  B  =  30. 


Termination  by  SLRT  with  A  = 
200  trials  per  parameter  pair 
* 

significant  at  p  <  0.05 


• 

. 
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TABLE  7 

CHARACTERISTICS  OF  RULE  R  WITH  AND  WITHOUT  DELAY 

u  •  ^ 


>  -1  _ 
2 

200 

ASN 

ITN 

OC 

x  -1 

1 

P 

Delayed 

Not 

Delayed 

Delayed 

Not 

Delayed 

Delayed 

Not 

Delayed 

600 

3 

152 

218 

33 

35 

1.000 

.998 

500 

2.5 

175 

231 

36* 

43 

.995 

.994 

400 

2 

212 

249 

47* 

58 

.915 

.949 

300 

1.5 

265 

273 

78* 

86 

.555 

.561 

266.67 

1.333 

248 

256 

77 

90 

.345 

.292 

200 

1 

210 

211 

- 

- 

.040 

.051 

150 

.75 

221 

225 

77 

81 

.255 

.284 

133.33 

.667 

219 

226 

70 

74 

.530 

.551 

100 

.5 

158 

178 

44 

46 

.920 

.943 

80 

.4 

124 

150 

26* 

33 

1,000 

.995 

66.67 

.333 

101 

134 

22 

26 

1.000 

.999 

x2_1  - 

10 

30 

3 

57 

56 

17 

17 

1.000 

.999 

25 

2.5 

64 

63 

21 

21 

1.000 

.993 

20 

2 

75 

79 

26* 

28 

.935 

.953 

15 

1.5 

103 

106 

40 

42 

.525 

.559 

13.333 

1.333 

98 

102 

41 

42 

.300 

.303 

10 

1 

77 

79 

- 

- 

.035 

.048 

7.5 

.75 

90 

98 

38 

41 

.295 

.297 

6.667 

.667 

94 

101 

38 

41 

.550 

.537 

5 

.5 

66 

73 

25 

27 

.940 

.948 

4 

.4 

47 

57 

17* 

21 

1.000 

.993 

3.333 

.333 

41 

53 

14 

20 

1.000 

.999 

Termination  by  SLRT  with  A  =  0.1,  B  -  30. 


200  trials  per  parameter  pair 
* 


significant  at  p  <  0.05 


' 
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APPENDIX 

COMPUTER  PROGRAMS 

I 

The  programs  used  to  derive  the  Monte  Carlo  results  of 
Chapter  IV  are  presented  below.  Program  MCI  simulates  the  R.,  rules, 
and  MC2  simulates  R^.  Both  programs  call  subroutine  DATA,  which 
allocates  storage  for  data  generated  by  subroutine  GGEXP  (from  the 
International  Mathematical  Subroutines  Library) . 

All  procedures  are  written  in  FORTRAN  IV,  for  execution  on 
the  University  of  Alberta  Computing  System. 

Following  are  the  input  variables  for  the  two  programs. 

1.  NRUNS  -  The  number  of  stopping  parameter  sets. 

2.  ISEED  -  The  "seed"  required  by  the  random-number  generation 

subroutines . 

3.  KRUNS  -  The  number  of  replications  desired  under  each  set  of 

parameters . 

4.  NALT  -  The  length  of  the  delay. 

5.  KARDS  -  The  number  of  pairs  of  exponential  parameters  to  be 

considered . 

6.  SL,  SLF,  SU,  SUF  -  A  stopping  parameter  set.  SL  corresponds  to 

A,  and  SU  to  B.  If  SLF  and  SUF  differ  from  SL  and  SU, 
respectively,  the  test  has  triangular  boundaries. 

7.  GAMMA  -  The  y  for  the  Ry  rules. 

8.  ST  -  The  maximum  number  of  trials  to  be  made.  This  was 

usually  taken  as  1200,  a  number  never  exceeded  in  the 
simulations . 

9.  C  -  The  critical  ratio  (p  ). 


non  non  non 
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C  PROGRAM  MCI 

C  SIMULATION  OF  3  (GAMMA)  WITH  TRUNCATION 

C 

COMMON  IS  FED,  X  (1  200,2)  ,  F  (2)  ,  MALL  (2)  ,  KLM  (2)  ,NP 
REAL  T  (2)  ,R  (2)  ,Y 
B5AL  X ,  E , EH  (2)  ,  AM  (2) 

INTEGER  KALL, KLM, N? 

INTEGER  D  (2)  ,N  (2)  ,  MT  (2) 

INTEGER  ISEED 

READ  (7,500)  NRUNS , KEUNS  , ISEED 
DO  7 5  KRUNCH= 1,HRUNS 

P.  E  A  D  ( 7 , 5  0  1 )  NALI,KAF.DS,SL,SLF,SU,SrJF,GAMMA,ST,C 
WRITE  (6,408)  SL,  SLF  ,  S'J  ,  SUF ,  ST ,  C  ,  N  ALT 
DO  75  L  = 1 , K ARDS 
READ  (7,  502)  E(1),E(2) 

WRITE  (6,407)  E  (1)  ,S  (2)  , GAMMA 
IN=2 

IF  (E  (1)  .LT.  S  (2)  )  I N=  1 
PTN=0 • 

A  S  N  =  0  • 

AM  (1)  =0. 

AM  (2)  =0. 

KH0  =  0 
KH  1  =0 
KH2=0 
A  2  -  0  • 

P2-0, 

M  A.XK-0 


INDIVIDUAL  RUNS  START  HERE 


DO  70  K  O  UN  T= 1 , K  R  U  N  S 
NN=0 

KALL  (1)  =1 
KALL  (2)  =  1 
KLM  (1)  =50 
KLM (2) =50 
DO  3  1=1,2 
NP=I 


SINGLE-RUN  INITIALISATION 

3  CALL  DATA 
DO  4  1=1,2 
T  (I) =0. 

D  (I)  -0 
N  (I)  =0 

4  ?  (I)  =1 . 

N?=  1 

CALL  GG'JB  (ISEED,  1  ,  Y) 

IF  (Y. LT*  0  • 5)  NP= 2 

SAMPLING  AREA 


1  N  (NP)  =N  (NP)  + 1 
NN=NN+  1 


■ 


nnn  nnnnn  onn 
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IF  (NN.GT. ST)  GO  TO  30 
IF  (N  (NP)  . LZ. KLK (NP) )  GO  TO  5 
CALL  DATA 
5  DO  15  J=1,2 


IF(N  (J) 

.  EQ. 

0) 

GO 

TO 

1  5 

K  =  N  (J) 

DO  10  I 

=  1,K 

II  (X  (I, 

J)  .E 

Q. 

o.) 

GO 

TO  10 

x  (I  /  J)  = 

X  (I, 

J) 

-1. 

IF(X  (I, 

J)  .G 

0.) 

GO 

TO  9 

D  (J)  =D  (J)  +1 

T  (J)  =T  (J)  +X 

(I 

,J)  i 

n. 

X (I,  J) =0. 

GO  TO  10 
9  T  (J)  =T  (J) +  1 . 

10  CONTINUE 

IF  (  (T  (1)  +T  (2)  )  .ZQ.O.  )  GO  TO  15 

?  (J)  =D  ( 3- J)  *ALOG  (C)  +  (D  (1)  +D  (2)  )  *  ALOG  (  (T  (1 )  +  T  (2)  )  / 

%  (C*T  ( 3- J)  +T  (J)  )  ) 

E  (J)  =  SXP  (E  (J)  ) 

15  CONTINUE 

X  =  I A  3  S  (D  (1)  -D  (2)  ) 

GO  TO  16 
12  n?=3-?:p 
GO  TO  1 

TERMINATION  RULES  APPLIED 

16  SB=SU+ (SUF-SU)  *NN/ST 
S A=SL+ (SLF-SL) *NN/ST 
A  ?  =  A  X  A  X  1  (R  (1)  ,E  (2)  ) 

IF  (  (AP.LT.SA)  .OR.  (AR.GT.  SB)  )  GO  TO  30 
ALLOCATION  RULES  APPLIED 
E  (GAME A)  RULE 

IF  (MAX  0  (D(1)  /  D(2)  )  .LT,  N  ALT)  GO  TO  12 
NP=2 

IF  (K.GT.  (GAXKA*NN)  )  GO  TO  20 
IF (T (1) +  T  (2) . SQ.0. )  GO  TO  20 
IF  (  (D  (1)/T(1)  )  .LE.  (D(2)/T(2)))  NP=1 
GO  TO  1 

20  IF  (D  (2)  .  GT.  D  (1)  )  NP=1 
GO  TO  1 

OUTPUT  ARIA 

30  IF  (D  (1)  .  EQ.  0)  EH  (  1)  =  1.  23456 
I F  ( D  ( 1 )  •  N E *  0 )  EH  (  1)=T  (1)/D  (1) 

IF  (D  (2)  .  EQ.  0)  EH  (2)  =  1.  23456 
IF  (D  (2)  .NS.  0)  EH  (2)  =T  (2)/D  (2) 

XT  (1)  =  N  (1)  -D  (1) 

XT  (2)  =  N  ( 2)  — D  (2) 

X  A  X  N  =  X A  X  0 {X AXNf NN ) 

WRITS  (3,4  14)  XT(1)  ,XT(2)  ,N(1)  ,N(2)  ,D  (1)  ,D  (2)  ,T  (1)  ,T  (2)  , 


r ' 

' 

. 
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C 

C 

C 


S ZH  { 1 )  ,  EH  (2) 

IF  (A?.,  LT.  SA)  GO  10  40 
IF  (NS.  GT.  ST)  GO  TO  4  5 
NP  =  1 

IF  (3  (2). GT.  ?.(*!))  NP  =  2 
GO  TO  (201,  20  2)  ,NP 

201  KK1  =  KH1  +  1 
GO  TO  50 

202  KE2  =  KH2  +  1 
GO  TO  50 

40  WRITS  (3,401)  NS 
KHO  =  KHO  +  1 
GO  TO  65 

45  AH=ALOG(EH(1)/2H(2)) 

IF  (ABS  (AH)  . LT. ALOG  (C) ) 
NP=  1 


IF  (AH.LT. 0. 0)  NP=  2 
50  WRITE  (3 ,402)  N  N  ,  N  P 
65  PTN=PTN  +  N  (IN) 

ASN=A  SN  +  N N 

A 2= A 2+ (FLOAT  (NN) )  **2 
P  2  =  P  2  +  N  (IN)  --2 
AM  (1)  =MT  (1)  +  AN  (1) 

AM  (2)  =MT  (2)  +  AM  (2) 

70  CONTINUE 


MULTI- RUN  SUMMARY 


GO  TO  40 


/ 


RUNS  =  FLOAT  (KRUNS) 

AS N= A  SN/EUN S 
PTN=PT  N/RUNS 

A  2=  (A2-3UNS*  (ASN**2)  )  /  (RUNS-1 . ) 

P2= (P2-3UNS*  (PIN**2) )/  (RUNS- 1 . ) 

AM (1) =  AM  (1) /RUNS 

AM  (2)  =  A  M  (2)  /RUNS 

WRITE  (6, 405)  ASN , PTN , A 2 , P 2 , KAXN 

WRITE  (6,411)  AM  ( 1 )  ,AM(2) 

75  WRITE  (6,  406 )  RUNS , KHO , K H 1 , KH2 
STOP 
C 

C  5XX=INPUT  4  XX=0  UTP  UT 

C 

50  0  FORMAT  (12,14,19) 

50  1  FORMAT  (14,1  2,  6? 6.  3,  F4.  1) 

502  FORMAT  (2F6.  2) 

407  FORMAT  ('1  MEAN  LIFE  ONE  =  ',F6.2,'  MEAN  LIFE  TWO  =  1 ,F6.2, 
S  2  0  X , 'GAMMA  =  »,F4.2) 

40  1  FORMAT  (*  AFTER  ',13,'  TFT  AT  NO  DIFFERENCE  FOUND') 

402  FORMAT  ('  AFTER  ',13,'  TRIADS  TREATMENT  ',11,'  CHOSEN') 

405  FORMAT  ('0  ASN  =  •  ,  F  5 .  1  ,  '  ITN  =  '  ,  F  5 .  1  ,  1  0  X  ,  '  V  AP.S  N  =  *,F10.2, 
S'  VARIN  =  ' ,  F10  •  2,  '  MAXN  =  ',14) 

40  6  FORMAT  ('  RUNS  =  ',F4.0,'  HO  ',14,'  HI  ',14,'  H2  ',24) 

408  FORMAT (' 1 ',////,' LOWER  TERMINATION  AT  ' ,F4.3,'  TO  ',F4.3, 


■ 


41  1 
4  1  4 


^///,  '0??rE  TERMINATION  AT  ,,F4.  1,'  TO  *  F4.1.///# 

5 •TRUNCATION  AT  1 , F6 . 1 ,//////, » CEITI C AL  RATIO  ,,?4.1, 
&///,  1  ALTERNATE  UNTIL  »,I4) 

FORMAT  ( 1  AVERAGE  NUMBER  OF  ACTIVES:  • , F6.  1 , 6X, F6.  1) 
FORMAT  (IX, 616,431 4.  5) 

END 


, 


no  o  non 
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C  PROGRAM  MC2 

C  SIMULATION  OF  R( E )  WITH  TRUNCATION 

C 

COMMON  ISEED,XC  1200,2  ),E(  2  ),  KALLC  2  ),KLM(  2  ),NP 
REAL  T(  2  )  ,  R(  2  ),  Y  ,  CH(  2  ) 

REAL  X  |  E  ,  E  E(  2  )  ,  A  M(  2  ) 

INTEGER  XA  LL , KLM , NP 
INTEGER  D(2),NC2),  MTC  2  ) 

INTEGER  ISEED 

REA D(  7, 500  )  NRUNS, KRUNS,  I  SEED 
DO  75  KRUNCH=1 , NRUNS 

RFAD(  7 , 501  )  NALT , KASDS , SL , SLF ,  SU, SUF, ST, C 
WRITE( 6,408)  SL , SLF , SU , SUF , ST, C , NALT 
DO  75  L= 1 , KARDS 
READ!  7, 502  )  EC  1  )  ,E(  2  ) 

WR  I  TEC  6,407)  EC  1  ),E(2) 

I  N=  2 

IFC  EC  1  ).  LT.EC  2  )  )  IN=1 
PTN=0. 

A  S  N  =  0  • 

AM(  1  )=9. 

A  *[(  2  )  =  0  • 

KH0  =  0 
KH  1=0 
KH2  =  9 
A2=  0  • 

P2  =  0  • 

MAXN=0 


INDIVIDUAL  RUNS  START  HERE 

DO  70  KOUNT=l , KRUNS 
NN  =  0 

KALLC  1  )=1 
KALLC  2  )  =  1 
KLMC  1  )  =  50 
KLUC  2  )  =  50 
DO  3  1=1,2 
NP=  I 


SINGLE-RUN  INITIALIZATION 


C 

C 

C 


3  CALL  DATA 
DO  4  1=1,2 
TC  I  )  =  0  • 

DC  I  )  =  0 
N(  T  )  =  0 

4  PC  I  )  =  1  • 

NP=  1 

CALL  GGUBC  ISEED,  l  ,  Y  ) 

IFC  Y. LT. 0. 5  )  NP  =  2 

SAMPLING  AREA 

1  NC  NP  )  =  N(  NP  )+  1 
NN=NN+l 


' 
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c 

c 


c 


c 

c 

c 

c 

c 


c 

n 

c 


IF<  NN.  GT.  ST  )  GO  TO  30 
IF(  N7(  NP  )•  LE  •  KLMl  NP  )  )  GO  TO  5 
CALL  DATA 
5  DO  13  J=  1 , 2 

IF(  N(  J  )•  SO.  0  )  GO  TO  15 
K=N(  J  ) 

DO  10  I  =  1,K 

IF(  XC  I  f  J  ).  EQ.O.  )  GO  TO  10 
XC  I ,  J  )  =  X(  I  ,  J  )- 1  • 

IFC  XC  I,  J  ).GE.O.  )  GO  TO  9 

D(  J  )  =  D(  J  )+  1 

T(  J  )  =  T(  J  )+XC  I,  J  )+l  . 

X(  T  f  J  )  =  0  • 

GO  TO  10 
9  T(J)  =  T(J)+1. 

10  CONTINUE 

IF(  (  T(  1  )  +  TC  2  )  ). EQ.O.  )  GO  TO  15 

R(  J  )  =  D(  3-J  )*ALOGC  C  )  +  (  D(  1  )  +  D(  2  )  )*ALOG(  (  T(  1  )  +  T(  2  )  )/ 
%C  C*T(  3-J  )  +  T(  J  )  )  ) 

P(  J  )  =  EXP(  R(  J  )  ) 

15  CONTINUE 

M=  I  ABS(  D(  1  )  —  D(  2  )  ) 

GO  TO  16 
12  NP=  3— NP 
GO  TO  1 

TERMINATION  RULES  APPLIED 

16  SB=SU+(  SUF-SU  )*NN/ST 
SA=SL+C  SLF-SL  )*NN/  ST 
AR  =  AMA X 1  (  R(  1  ),  R(  2  )  ) 

IF(  (  AR.LT.SA  ).OR.(  AR.GT.SB  )  )  GO  TO  30 
ALLOCATION  RULES  APPLIED 
RC  E  )  RULE 

IFC  MAX0C  DC  1  ),D(  2  )  ).  LT.NALT  )  GO  TO  12 
MD=MIN0(  D(  1  )  ,D(  2  )  > 

I  F(  MD.  GT.  0  )  TO  TO  19 
TF(  M  D •  E Q  •  D (  1  )  )  NP=1 
IFC  MD. EO.D(  2  )  )  NP=2 
GO  TO  1 

1  9  B=(  T(  1  )*D(  2  )  )/ Tl  1  )*D(  2  )  +  T(  2  J(  1  )  ) 

CH  (  1  )  =  A  B  S  (  B  — (  N(  1  )+l  .  0  )/(  N(  1)+M(2)+1.0  )  ) 

CH(  2  )  =  A B S (  E— N(  1  )/(  N(  1  )  +  N(  2  )+ 1 . 0  )  ) 

DO  20  1=1,2 

20  IF(  AM  INK  CH(  1  ),CF(  2)  ).EQ.CH(  I  )  )  NP=I 
GO  TO  1 

OUTPUT  AREA 

30  I  F(  D(  1  )  •  EQ  •  0  )  Eh'(  1  )=1 . 23456 
I  F (  D(  1  )  •  NE  •  0  )  E H<  1  )  =  T(  1  >/D(  1  ) 

IFC  D(  2  ).  EQ.O  )  EH(  2  )=1 . 23456 
IFC  DC  2  )  .NE.O  )  EHC  2  )  =  TC  2  )/DC  2  ) 
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MT(  1  }  =  N(  1  )  — Dl  1  ) 

MT.  2  >  =  N(  2 )  — Dt 2  ) 

M A  XN  =  M A X  0 (  M A XN , X  N  ) 

WRITE(  3,414)  MTl  1  )  ,  MTt  2),N(1  ),N(  2  ),  D(  1  ),  D(  2  ),  T(  1  ) ,  T(  2  )  , 
£FH(  1  )  ,  EH(  2  ) 

IFtAR.LT.SA)  GO  TO  40 
TFtNN.GT.ST)  GO  TO  45 
N?=l 

I  Ft  R  t  2  )  •  GT  •  R (  1  )  )  NP=2 
GO  TO  (  201 ,202 ), NP 

201  KHl  =  KH1  +  1 
GO  TO  50 

202  KH2  =  KH2  +  1 
GO  TO  50 

4  0  WRITEt  3,401)  NNT 
KH0  =  XKO  +  1 
GO  TO  65 

45  AH=ALOG(  EH(  1  )/EH(  2  )  ) 

IFt  ABSt  AH  ).  LT.  ALOGt  C  )  )  GO  TO  40 
NP=1 

I F( AH. LT. 0. 0  )  NP  =  2 
50  WRITEt 3,402)  NN, NP 
65  PTN  =  PTN+Nl  IN  ) 

ASN  =  ASN+NN 

A2  =  A2+t  FLCATt NN )  )**2 
P2=P2+Nt  IN  )**2 
AM(  1  )  =  MTt  1  )+AM(  1  ) 

AMt  2  >  =  MTt  2  )  +  AMt  2  ) 

70  CONTINUE 
C 

C  MULTI -RUN  SUMMARY 

C 

RUNS  =  FLOAT(  KRUNS  ) 

ASN= ASN/ RUNS 
PTN=  PTN/ RUNS 

A2  =  (  A2-RUNS*t  ASN**2  )  )/t  RUNS-1  .  ) 

P2=(  P2  -  RUN  S^(  PTN**2  )  )/(  RUNS-1  .  ) 

AM(  1  )  =  AM(  1  )/  SUNS 
AM(  2  )  =  AM(  2  )/ RUNS 

WR ITE(  6, 405  )  ASN , PTN , A2 , P2 , MAXN 
WRITEt  6, 411)  AMt  1  )  ,  AMt  2  ) 

T5  WRITEt 6,406)  RUNS , KHO , KHl , KH2 
STOP 
C 

C  5XX= I NP  UT  4 XX  =  OUTPUT 


500 

FORMATt 

12, 14,19) 

501 

FORM \Tt 

14, 12  , 6 F6 

.3,  F4 

.  1  ) 

50  2 

FORMATt  2F6 . 2  ) 

407 

FOPMATt 

« 1  MEAN 

LI  FE 

ONE  =  '  , 

F  6 • 2 ,  *  MEAN  LIFE  TWO  =  ',F6.2) 

40  1 

FORMATt 

'  AFTER  * 

,13,  • 

TRIALS 

NO  DIFFERENCE 

FOUND'  ) 

402 

FORMATt 

•  AFTER  ' 

,13,  ' 

TRIALS 

TREATMENT  ',11 

,'  CHOSEN*) 

405 

FORMATt 

*0  ASN  = 

*  ,  F5. 

1  ,  *  I TN 

=  '  ,F3. 1 , 1  OX,  * 

VARSN  =  *  ,  F 1 0 • 2 , 

■ 

. 
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P '  VAR IN  =  •,FlC.2f«  MAXN  =  »  ,14) 

-JO  A  FORMAT!'  RUNS  =  »,F4.0,'  HO  «f  I4f  '  Hi  ',I4,»  H2  ',14) 

4 0 S  FORM AT(  '1', ////,' LOWER  TERMINATION  AT  »,F4.3,'  TO  »,F4.3, 

P/  /  / ,  • UPP  ER  TERMINATION  AT  ',F4.1,'  TO  ',F4,l,///t 
S'  TRUNCATION  AT  ' fF6« 1 ,//////,• CRITICAL  RATIO  ',F4.1, 

P/// ,  ' ALTERNATE  UNTIL  ',14) 

411  FORMAT!  '  AVERAGE  NUMBER  OF  ACTIVES:  '  , F 6 . 1 , 6X , F6 . 1  ) 

414  FORMAT!  1 X , 6  I  6 , 4 E 1 4 . 6  ) 

END 


■ 

. 
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SUBROUTINE  DATA 

C  SUBROUTINE  FOR  GENERATION  AND  ALLOCATION  OF  DATA 

COMMON  I  SEED  *  Y(  1200,2  )  ,  E(  2  ),KALL(  2  ),KLM(  2  )  ,  NP 
REAL  Y,E,YV( 50) 

INTEGER  KALL, KLM ,NMP, I  SEED 
CALL  GGEXP(  I3EED , E( NP  ), 50, YV) 

NNN=F0*(  KALL(  NP  )-l  )+l 
NNNN=NNN+49 
DO  1000  I  =  N N  N , N  N  N  N 
J=MOD(  1,50) 

IF(J.EQ.O)  J  =  50 

1  000  Y(  I , NP  )  =  YV { J  ) 

KLM(  NP  )=50*KALL(  NP  ) 

KALL(  NP  )  =  KALL(  NP  )+l 

RETURN 

END 


' 
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